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Abstract 

This paper describes the measurement of elliptic flow of charged particles in 
lead- lead collisions at = 2.76 TeV using the ATLAS detector at the 

Large Hadron Collider (LHC). The results are based on an integrated lumi- 
nosity of approximately 7 /xb~^. Elliptic flow is measured over a wide region 
in pseudorapidity, < 2.5, and over a broad range in transverse momentum, 
0.5 < pt < 20 GeV. The elliptic flow parameter V2 is obtained by correlating 
individual tracks with the event plane measured using energy deposited in the 
forward calorimeters. As a function of transverse momentum, V2{pt) reaches a 
maximum at pt of about 3 GeV, then decreases and becomes weakly dependent 
on pt above 7-8 GeV. Over the measured pseiidorapidity region, V2 is found to 
be only weakly dependent on r], with less variation than observed at lower beam 
energies. The results are discussed in the context of previous measurements at 
lower collision energies, as well as recent results from the LHC. 
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1. Introduction 

The measurement of collective phenomena in nuclear collisions at high en- 
ergies has been a subject of intensive theoretical and experimental studies. 
Anisotropic flow, which manifests itself as a large anisotropy in the event-by- 
event azimuthal angle distribution of produced particles, is generally understood 
to be a consequence of the spatial anisotropy of the initial energy deposition from 
nucleon-nucleon collisions in the overlap of the colliding nuclei. Anisotropics in 
the initial energy density are converted into final state momentum anisotropics 
via strong rescattcring processes which induce pressure gradients, following the 
laws of relativistic hydrodynamics. Consequently, azimuthal anisotropics are 
sensitive to the initial state and its subsequent dynamical evolution. 
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Anisotropic flow is commonly studied by measuring the Fourier coefficients 
(vn) of the azimuthal angle distributions of the emitted particles. The second 
harmonic, V2, referred to as "elliptic flow", is the most extensively studied as 
it most directly relates the anisotropic shape of the overlap of the colliding 
nuclei to a corresponding anisotropy of the outgoing momentum distribution 
(for a review, see Ref. pLj). Elliptic flow has been measured over a wide range of 
energies, collision systems, and collision centralities by all of the RHIC heavy ion 
experiments [2l|3l|4l[5] and several experiments at lower energies (see Ref. [1]). 
Predictions for V2 at the LHC energy varied widely, covering all possibilities 
from a strong rise, no change, or even a decrease of V2 [B] relative to lower 
energy collisions. Measurements of V2 for inclusive charged-particles from the 
ALICE experiment [7] indicate that, integrated over pt, it increases by about 
30% from RHIC to LHC energies. However, ALICE also observed that V2{pt) 
for inclusive charged particles was identical with RHIC results for the same 
collision centrality (or impact parameter) up to pj^ — A GeV. This implies that 
the observed rise is driven primarily by an increase in the average transverse 
momentum with the higher collision energy. 

In this letter, we present a measurement of the elliptic flow of charged parti- 
cles in lead- lead collisions at y/s^ = 2.76 TeV with the ATLAS detector at the 
LHC. The elliptic flow is measured in the pseudorapidity region I77I < 2.5 over 
the full azimuthal range < ip < 27r, for transverse momentgj^O.S < < 20 
GeV. This allows stringent tests of the applicability of hydrodynamics in the 
LHC energy regime, and provides information on the transition between low pt, 
where hydrodynamics is expected to dominate, and higher px, where particle 
production is expected to stem from the fragmentation of jets modifled by the 
hot, dense medium ,8,. 

2. ATLAS detector and trigger 

The ATLAS detector [S] is well suited for measurements of azimuthal an- 
isotropies over a large pseudorapidity range. The relevant detectors for this 
analysis are the inner detector (ID) and forward calorimeter (FCal). The ID is 
contained within the 2 T field of a superconducting solenoid magnet, and mea- 
sures the trajectories of charged particles in the pseudorapidity region \r]\ < 2.5 
and over the full azimuthal range. The precision silicon tracking detectors con- 
sist of pixel detectors (Pixel) and a semiconductor microstrip tracker (SCT). 
In the "barrel" region, these are arranged on cylindrical layers surrounding the 
beam pipe, while in the "endcap" regions they are mounted on disks perpendic- 
ular to the beam axis. A charged particle typically traverses three layers of the 



^ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction 
point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the y axis points upward. CyUndrical 
coordinates (r, 0) are used in the transverse plane, being the azimuthal angle around the 
beam pipe. The pseudorapidity is defined in terms of the polar angle 8 as rj = — lntan(9/2). 
Transverse momentum and energy are defined as pT = psinS and i?x = Esin9, respectively. 
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Pixel detector and four double-sided layers of the SCT. The silicon detectors 
are surrounded by a transition radiation tracker (TRT), composed of drift tubes 
and covering up to \ri\ = 2. 

The FCal covers a pseudorapidity range 3.2 < |ry| < 4.9. It uses tungsten and 
copper absorbers with liquid argon as the active medium, with a total thickness 
of about 10 interaction lengths. This analysis uses the energy deposition in 
the entire FCal for the centrality determination, while for the reaction plane 
measurement only the energy deposition in the first sampling layer of the FCal 
(Layer 1) is used, as doing this was found to minimize the effect of fluctuations 
on the reaction plane measurement. 

The trigger system has three stages, the first of which (Level-1) is hardware- 
based, while the later stages (Level-2 and Event Filter [3]) are based on software 
algorithms. The minimum-bias Level-1 trigger used for this analysis requires sig- 
nals in either the two sets of minimum-bias trigger scintillator (MBTS) counters, 
covering 2.1 < < 3.9 on each side of the experiment, or the two zero-degree 
calorimeters (ZDC), each positioned at \z\ = 140 m relative to the centre of 
the detector, detecting neutrons and photons with \r]\ > 8.3. The ZDC Level-1 
trigger thresholds were set just below the single neutron peak on each side. The 
MBTS trigger was configured to require at least one hit above threshold from 
each side of the detector. A Level-2 timing requirement on signals from the 
MBTS was then imposed to remove beam backgrounds, while the ZDC had no 
further requirements beyond the Level-1 decision. The Event Filter was not 
needed for the minimum-bias triggering and was run in pass-through mode. 

3. Event selection and reconstruction 

The lead-lead data set analysed here corresponds to an integrated luminos- 
ity of approximately Lint = 7 ^h~^. Three main event selection requirements 
were applied offline to reject both non-collision backgrounds and Coulomb pro- 
cesses, in particular highly-inelastic photonuclear events. First, an offline event 
selection required a time difference |At| < 3 ns between the positive and neg- 
ative 77 MBTS counters as well as a reconstructed vertex in order to suppress 
non-collision backgrounds. Second, a coincidence of the ZDCs at forward and 
backward pseudorapidities was required in order to reject a variety of back- 
ground processes, while maintaining high efficiency for non-Coulomb processes. 
Finally, in this analysis only events with a vertex with |zvtx| < 10 cm were used. 
Simulations show the vertex algorithm to be essentially 100% efficient for the 
event sample considered here. Pile-up events, defined as additional minimum 
bias events in the same bunch crossing, are expected to be present at the 10"'' 
level and so are negligible. In total, approximately 4 x 10'' events passed the 
selection criteria. 

Tracks were reconstructed within the full acceptance of the inner detector. 
To improve the reliability of the ID track reconstruction in the tracking environ- 
ment in heavy ion collisions, the track quality requirements are more stringent 
than those defined for proton-proton collisions [10] . Tracks are required to have 
at least eight hits in the SCT, at least two Pixel hits and a hit in the Pixel layer 
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closest to the interaction point. A track must have no missing Pixel hits and 
at most one missing SCT hit, where such hits are expected. Finally, the trans- 
verse and longitudinal impact parameters with respect to the vertex (|do| and 
1^0 sin 0|) were each required to be less than 1 mm. These additional require- 
ments were made to improve the purity of the track sample. The inefficiency of 
this selection is driven by the loss due to hadronic interactions in the detector 
material, which increases with [in]- Thisresuhs in an additional inefficiency 
of approximately 15% at \r]\ > 1 compared to the central region of the detector. 

However, the results shown here are found to be insensitive to the absolute 
tracking efficiency (discussed below) , and the effect of the efficiency decrease at 
high 1 77 1 is minimized when measurements are performed in limited transverse 
momentum and pseudorapidity intervals. 

The tracking performance has been studied in detail by comparing data 
to Monte Carlo simulations based on the HIJING event generator [TT] and a 
full GEANT4 [T^ simulation of the detector |T3j. In general the simulated 
distributions of the number of Pixel and SCT hits on tracks describe the data 
well, particularly after reweighting the simulated momentum distribution to 
account for the differences in the charged particle spectrum reconstructed in 
data and HIJING. Monte Carlo calculations show that the tracking efficiency 
for charged hadrons in this analysis is about 72% near = in central collisions, 
lower than in proton-proton collisions due to the more stringent requirements 
and the higher occupancy in the SCT. Fake tracks from random combinations of 
hits are generally negligible, e.g. reaching only 0.1% in Irj] < 0.3 for the highest 
multiplicity collisions, although the rate of fake tracks increases slightly with 
increasing \r]\. 

4. Data analysis 

In order to systematically select various geometries of the initial state, the 
data were analysed in centrality intervals defined by selections on FCal T,Et, the 
total transverse energy deposited in the FCal (always stated at the electromag- 
netic energy scale, which does not correct for the response of the calorimeter to 
hadrons). These intervals are expressed in percentiles of the total inelastic non- 
Coulomb lead-lead cross section (0-10%, 10-20%,..., 70-80%) with the most 
central interval (0-10%) corresponding to the 10% of events with the largest 
FCal Ei?T- The measured FCal Y,Et distribution for a subset of the data (with 
Lint approximately 200 mb^^), taken with a less restrictive primary trigger than 
used for the bulk of the data and used for the calibration procedure described 
below, is shown divided into centrality intervals in Fig. [T] 

To establish the fraction / of the total non-Coulomb inelastic cross section 
selected by our trigger and event selection criteria, we have performed a convo- 
lution of FCal "EEt distributions measured in proton-proton data at ^/s = 2.76 
TeV with a full Monte Carlo Glauber calculation [14] . The calculation assumes 
the number of effective proton-proton collisions per lead-lead event, A'", scales 
according to the "two-component model" (from e.g. Ref [H]). This model 
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Figure 1: Measured FCal SBy distribution divided into 10% centrality intervals (black). 
Proton-proton data at y/s = 2.76 TeV, convolved with a Glauber Monte Carlo calculation 
with X = 0.088 (grey), as described in the text. 



combines the number of participants ( A'part i the number of nucleons which in- 
teract inelastically at least once) and the number of binary coUisions (A'coii) as 
= (1 — x)^^^^ -{-xNcoW- In this approach, the only free parameter is x, which 
controls the relative contribution of A'part and NcoW- The best description of the 
data is found to be for x = 0.088. The value of / and its uncertainty is estimated 
by systematically varying the effect of trigger and event selection inefficiencies 
as well as backgrounds in the most peripheral FCal TjEt interval to achieve the 
best agreement between the measured and simulated distributions. Using this 
analysis of the FCal YiEt distribution, the fraction of the total cross section 
sampled by the trigger and event selection has been estimated to be 98%, with 
an uncertainty of 2%. This is similar to estimates given in a previous ATLAS 
publication [TB]. The FCal Y^Et ranges defined from this subsample have been 
found to be stable for the full data set, both by counting the number of events 
and by measuring the average number of reconstructed tracks in each interval. 
The 20% of events with the smallest FCal YjEt are not included in this analysis, 
due to the relatively large uncertainties in determining the appropriate selection 
criteria. 

The final state momentum anisotropy can be quantified by studying the 
Fourier decomposition of the azimuthal angle distribution |17) : 

^d^N _ 1 (fN _ 1 E cPN Lj^^S^v cos[7i(0 ^ )]] (1) 
dp^ Pt d(f>dpT^dy 27rpT P dp^dij \ " n i , 

where ?/, pt and </> are the rapidity, transverse momentum, and azimuthal angle 
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of final-state charged particle tracks and denotes the azimuthal angle of the 
n-th order reaction plane. In more peripheral events, ^2 is close to ^rp, the 
reaction plane angle, defined by the impact parameter (b, the vector separation 
of the barycentres of the two nuclei) and the beam axis {z). In more central 
events, ^'2 primarily reflects fluctuations in the initial-state conflgurations of 
colliding nucleons. This analysis was confined to the second Fourier coefficient 
{n = 2), V2 = (cos [2(0 — <i>iip)]), where angular brackets denote an average 
first over particles within each event relative to the eventwise reaction plane, 
and then over events. 

In this analysis, the n = 2 event plane is determined from the data on an 
event-by-event basis, according to the scheme outlined in Ref. [17^: 

^'^2'''' l,Ei^^T'-.cos(20,)j' 

where sums run over tower transverse energies E^™'^'^ as measured in the first 
sampling layer of the forward calorimeters, with each tower covering Arj x A(j) = 
0.1 X 0.1. The tower weights, Wi = Wi{(t)i,rii), are used to correct for local 
variations in detector response. They are calculated in narrow A?7 slices (A77 = 
0.1) over the full FCal range in such a way as to remove structures in the 
uncorrected (/) distributions of in every A?7 slice. The final results of 

this analysis are found to be insensitive to the weighting, and results obtained 
with all Wi = 1 were consistent with those reported here, and well within the 
systematic uncertainties estimated below. 

The correlation of individual track azimuthal angles with the estimated event 
plane is shown in Fig. [2] for tracks with px = 1 — 2 GeV. There is a clear sinu- 
soidal modulation at all centralities. The modulation is largest in the 20-50% 
centrality intervals, and decreases for the more central and peripheral events. In 
the centrality intervals where the correlation is strongest, the correlation does 
not follow a perfect l-\-a cos(20) form, indicating significant contributions from 
higher order harmonics. However, in this letter we rely on the orthogonality of 
the Fourier expansion and do not extract the other coefficients. To verify that 
this does not bias the measurement, we have extracted V2 from a fit containing 
all Fourier components w„ up to n = 6, and found V2 values consistent with the 
results extracted below. The odd amplitudes are found to be consistent with 
zero, as expected when measuring odd harmonic functions relative to ^2 |17j . 

The measured values of V2 are generally underestimated because of the finite 
experimental resolution in extracting the event plane angle. The event plane 
resolution correction factor, i?, was obtained using the subevent technique, also 
described in Ref. [T7]. Two "subevents" are defined in each event, one each 
in the forward and backward 77 directions. For the measurement of the event 
plane using the FCal, the first sampling layer on the positive 77 side was selected 
as subevent "P", with a corresponding subevent 'W formed for negative 77. 
The resolution correction for the event plane measured by each subevent was 
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calculated as a function of FCal Ei?T according to the formula 



i?(I]i?T) = v'(cos[2(M^r-*2^)]), (3) 

where angular brackets denote an average over all events in a FCal 'SEt interval. 
The left panel of Fig. [s] shows the distribution of the difference — ^t^. The 
right panel shows the FCal Y,Et dependence of the resolution correction for 
the event plane determined using the full FCal Layer 1 as well as a reduced- 
acceptance version used in the systematic studies discussed below. 

The final, resolution-corrected, V2 is calculated in intervals of centrality, 77 
and pt as 

-. events tracks 

V2{v,Pt)-^Y. .pii^ E ^.cos[2(<^.-<f)], (4) 

where iV^'ot denotes the total number of the reconstructed tracks in a given 
centrality, rj and px range, and the Ci are weights similar to the Wi for tracks, 
designed to flatten the (p distribution in a small A77 slice. For '^^j^ (the event 
plane for event j) we take the event plane measured in the opposite 77 hemisphere 
(i.e. "P" at positive 77, or 'W at negative 77) to each track with azimuthal angle 
Using the track in the opposite hemisphere maximizes the pseudorapidity 
gap between the reaction plane estimate and the V2 estimate (IA77I > 3.2), 
minimizing potential non-flow correlations between them. 

The systematic uncertainty on V2 as a function of V E^nd centrality was 
evaluated by varying several different aspects of the analysis procedure. 

• The resolution correction was changed by limiting the FCal acceptance to 
a smaller range in pseudorapidity. 

• Tighter tracking requirements were applied (both |c?o| and |zosin0| less 
than 0.5 mm, instead of the nominal 1 mm requirement). 

• Results were compared using negatively and positively charged tracks. 

• Results were compared between V2 measured at positive and negative pseu- 
dorapidities. 

• Results were studied as a function of time during the heavy ion run. 

Additional sources of systematic uncertainties were examined, including the 
following: Deviations from zero of (sin(2[(/) — \E'2])), which are sensitive to resid- 
ual biases in the reaction plane determination and detector non-uniformities, 
were measured. Monte Carlo studies were performed based on HIJING, with 
a special procedure applied to the generated particle azimuthal angles so as 
to simulate elliptic flow (from Ref. [H]), with a magnitude extrapolated from 
RHIC data. Deviations from the flow induced at the generator level were ob- 
tained by applying the same analysis procedure to the simulated data as with 
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Centrality 


0-10% 


40-50% 


PT [GeV] 


0.8-0.9 


2.4-2.7 


8-10 


0.8-0.9 


2.4-2.7 


8-10 


Smaller r] acceptance of event 


0.6 


1.2 


5.7 


0.7 


0.7 


2.0 


plane determination 














Residual deviation from zero of 


0.7 


0.6 


0.4 


0.5 


0.7 


1.2 


sine terms 














Varying tracking cuts 


0.4 


0.1 


1.7 


0.1 


< 0.1 


0.2 


Negative vs positive tracks 


0.5 


0.3 


3.3 


0.3 


0.1 


1.6 


Asymmetry with respect to rj 


0.1 


0.1 


0.2 


< 0.1 


< 0.1 


0.1 


reflection 














Time dependence 


0.2 


0.2 


Monte Carlo reconstruction 


1.2 


1.2 


1.2 


0.3 


0.3 


0.3 


Total systematic error 


1.6 


1.9 


6.9 


1.0 


1.0 


2.9 



Table 1: Principal systematic uncertainties (stated as a percentage of the value of V2) on the 
V2 measurement for three p-p intervals and two centrality intervals, all for |»?| < 1. 



real data. The event plane determined from the reconstructed tracks was also 
investigated as an independent cross-check on the FCal reaction plane. In this 
case, for the tracks with positive (negative) 77 the event plane determined in the 
negative (positive) rj subevent was used. The uncertainty in the fraction of the 
total inelastic cross section sampled by our trigger and event selection criteria 
gives an overall scale uncertainty on V2, ranging from 1% in central events up 
to 5% in peripheral events. 

Deviations in individual contributions from the baseline results have been 
quantified as relative systematic uncertainties on V2 (in percent), which are 
listed in Table [l] for several centrality and pt intervals, all for \ri\ < 1. The 
different components have been added in quadrature and expressed as la point- 
to-point systematic uncertainties. Note that the somewhat large increase in the 
scale of the uncertainties from moderate to high can be partly attributed to 
the limited track statistics at high px- It should also be pointed out that the 
systematic uncertainties only include those associated with the measurements 
themselves; no attempt is made to disentangle the potential contributions from 
non-flow effects, since their nature is not yet fully understood. 

5. Results 

The top panel of Fig. |4] shows the f 2 dependence on px for eight 10% cen- 
trality intervals and for tracks with I77I < 1. It is observed that all centrality 
intervals show a rapid rise in V2{pt) up to = 3 GeV, a decrease out to 7-8 
GeV, and then a weak px-dependence beyond 9-10 GeV. The same trends are 
also seen for 1 < It^I < 2 (Fig. [i] middle) and 2 < \ri\ < 2.5 (Fig. [i] bottom). 

The pseudorapidity dependence of V2 is shown in Fig. [5] The top row shows 
the centrality and ry dependence of V2{'i],p^) for five px intervals, which char- 
acterize the trend shown in Fig. [4j and the four most-central intervals. The 
bottom row shows the same information for the four most-peripheral intervals. 
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It is observed that V2 depends very weakly on rj over the measured pseudora- 
pidity region. In the two lowest intervals, below 1.2 GeV, V2 drops by about 
5-10% over the range jr/j = — 2.4. At higher transverse momenta, a decrease 
on the order of few percent can be seen, although, due to the large point-to- 
point errors, a flat 77 dependence cannot be excluded. This is in contrast to the 
strong variation in V2{vi) observed by the PHOBOS experiment at ^snn = 200 
GeV [IHj, which drops by approximately 30% between = and ry = 2.5. 

Fig. |6] shows W2(pt) for |ry| < 1 in the 40-50% centrality interval compared 
to data from the LHC (ALICE, from Ref. [7]) as weU as from RHIC (STAR [19] 
and PHENIX [20J) with a centre-of-mass energy a factor of nearly 14 lower. The 
ALICE and STAR data are shown for the second cumulant W2{2}, which gives 
results closest to the event-plane method used in this analysis. The PHENIX 
data are obtained with a similar method as ATLAS, but with V2 measured only 
for identified 7r° hadrons, detected through their two-photon decay mode. It is 
observed that all of the data sets are quite similar as a function of p^, both at 
lower pt (ALICE and STAR) and even at higher px, within the limited statisti- 
cal precision of the PHENIX data. The observation of similar V2 at low has 
been noted recently [7| , and has been reproduced using hydrodynamical simula- 
tions assuming the same shear viscosity to entropy density ratio but initialized 
at a higher energy density. However, the similarities at high px will require 
additional theoretical study to see if they are consistent with the differential 
energy loss of jets in the hot, dense medium. 

6. Conclusions 

Elliptic flow measurements in lead- lead collisions at — 2.76 TeV ob- 

tained with the ATLAS detector are presented for an integrated luminosity of 
approximately 7 /ib^^. These results represent the first measurement of V2 over 
a broad range in r/ and px at the LHC energy. As a function of transverse mo- 
mentum, at all |r/|, V2 rises rapidly up to — S GeV, decreases somewhat less 
rapidly out to px = 7-8 GeV, and then varies weakly out to 20 GeV. Over the 
measured pseudorapidity region, \ri\ < 2.5, V2 is found to be only weakly depen- 
dent on 77, with less variation than observed at lower beam energies. Comparison 
of the 40-50% interval with lower energy data shows little change both at low 
and high px- These results provide strong constraints on models which aim 
to describe the dynamical evolution of the system created in ultra-relativistic 
heavy ion collisions. 
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Figure 2: Distribution of the azimuthal angle of individual tracks relative to the measured 
event plane, in eight centrality intervals. These distributions are meant to illustrate the 
observed correlation relative to the event plane, and are not used in the quantitative estimates 
of V2- The curve is a fit to 1 -|- cos(n</>) up to n = 6. 
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Figure 3: (left) Distribution of the difference between the event planes at positive and 
negative r/ obtained using Layer 1 FCal towers, both with full and half acceptance, (right) 
FCal TiEx dependence of the resolution correction for event planes from Layer 1 FCal towers 
in full acceptance (full symbols) and half acceptance (open symbols). 
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Figure 4: Elliptic flow V2{pt) as a function of p-p for eight 10% ccntrality intervals, for 
Pt from 0.5 to 20 GeV, and for three ranges in pscudorapidity (|r;| < 1, 1 < |r;| < 2 and 
2 < < 2.5). Error bars show statistical and systematic uncertainties added in quadrature. 
The arrows indicate where the value of V2 does not fit within the chosen plot scale, due to 
large statistical fluctuations. 
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Figure 5: Pscudorapidity dependence of V2{pt- '/) for 0.5 < px < 20 GcV in five p-p intervals 
and 10% centrality intervals. Error bars show statistical and systematic uncertainties added 
in quadrature. 
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Figure 6: V2 vs. px at |?7| < 1 in the 40—50% centrality interval, compared to previous exper- 
imental data: ALICE t'2{2}[7] for inclusive charged particles, PHENIX I20| t>2 for identified 
tt", and STAR data on i>2{2} for inclusive charged particles for the 40-60% interval [19] . 
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